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1. Introduction 

A precise description of the neutrino-nucleus interaction is crucial for the study the neutrino 
properties such as CP violating phase and mass hierarchy [||, [|, [|,|4j] when one extracts the neutrino 
flux from the neutrino data. The quasi-elastic (QE) reaction is the key mechanism to reconstruct the 
neutrino flux from the data of neutrino-nucleus reaction. Neutrino energy (E v ) can be determined 
from the muon scattering angle(0) and cncrgy(E u ) by using the kinematics of the reaction of free 
nucleon at rest. 


EfjM N — m^/2 


(1.1) 


Mn — E^ + p^ cos 0 


Since the reaction takes place inside nuclear many body system, various initial and final state in¬ 
teractions can spoil the useful relation of Eq. O- The observed quasi-elastic like events are 
mixture of the various reaction mechanism, which contribute the reaction amplitudes and, in prin¬ 
ciple, interference among various mechanism exists. The central issue of the current theoretical 
research of neutrino-nucleus reaction is to improve the accuracy of the estimation of various nu¬ 
clear effects^, []| by testing the model from the neutrino data. 

Because of the rather wide energy band of the neutrino flux in the current long baseline ex¬ 
periments, the observed cross section of neutrino reaction (< a >) is a result of neutrino flux(<t> v ) 
average of the cross section (a). 



( 1 . 2 ) 


One can regard the extraction of neutrino flux <t> v from the data < a > as an inverse problem 
to obtain <t> v from < a > provided the cross section a is well under control. In this report we 
will examine a method to solve the inverse problem, where the neutrino flux can be extracted 
without using the formula ( |Tt| ) nor the separation of ’true’ QE events. For this purpose, we use the 
maximum entropy method (MEM) which are widely used in the field of condensed matter physics 
and Lattice QCD[|7|, [S[ |hj|. 

We briefly outline the application of MEM to extract neutrino flux from data in section 2. To 
examine the usefulness of the approach, we generate pseudo data of neutrino reaction. The pseudo 
data are generated by using a model of neutrino reaction and neutrino flux, which are explained in 
section 3. Here, as an example of the pseudo data, we use double differential cross section of muon 
for the inclusive neutrino nucleus reaction. Results on how well the neutrino flux can be extracted 
from the pseudo ’data’ without assuming QE mechanism is shown in section 4. 

2. Maximum Entropy Method 

In this section we briefly explain the maximum entropy method. The most plausible neutrino 
flux from the experimental data based on the Bayes’ theorem is given by the functional 

integral of the neutrino flux 


qMEM 


J [tl(pft>I J \<P\D.I\ 


( 2 . 1 ) 
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where /’[<f>|D,7] is a conditional probability of neutrino flux <f> for given data of neutrino reaction 
D = {< a > eX p) and the prior information on the neutrino flux I. I is called as default model of 
neutrino flux. Introducing auxiliary variable a, we can write the above formula as 

^MEM = I da <pap{ a | 5 ) /] ) ( 2 . 2 ) 

where d>“ is neutrino flux for given a, data D and prior information / given as, 

= j[d<t>]<i>P[<i>\a,DJ] °= j[dd>]<i>P[D\<i>,aJ]P[d>\a,I]. (2.3) 

Here P[D |d>, a,I] and Pdfcx.I] are called as likefood function and prior probability, respectively. 
The likefood function is written as 

P[D\^,I\ = i-exp^-i* 2 ), (2.4) 

with 

^2 _ ^ (A - A)~ 

i=i °F 

Z x = Y\J^f. (2.5) 

i=i 

Here D\ and D\ are cross sections < a > exp and < a >. 

The prior probability can be expressed as 

P[<E>|a,/] = ^-exp(aS), (2.6) 

Zs 

where S is the Shannon-Jaynes entropy, 

S = (2.7) 


mi is default model of neutrino flux at neutrino energy £). Combining the likefood function and 
the prior probability, neutrino flux <l>“ is the maximum probability of P[<t>\a,D,I] °c with 

Q{&) = &S-\l 2 - 

Finally the a dependent flux is integrated with the probability P[a\D,I], which has a sharp 
peak as a function of a written as 


P[a\D,I] 


/] 


P[a\I\j [d< I>] 


1 c 6(0) 
Z S Z X 


( 2 . 8 ) 


Assuming constant P[a\I], we obtain the neutrino flux (p^ EM by integrating a around the shark 
peak of P[a\D,I\ in Eq. |2.2| . 

One can also estimate the errors of the extracted neutrino flux within MEM. For given 

a, the average flux < <f>« >s in energy region S of neutrino is defined as 


{(p a )s = J S dEv^(E v ) 


J s dE v 


(2.9) 
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The covariance of (<f>“) 5 can be written as 

;(<5<I>“) 2 ) s = ^ dE v dE' v 8<f>(E v ) < f>(E' v )P[^>\D,a,I]/ ^dE v dE' v 

8 2 Q 


Is 2 

— j dE v dE{, 


s 2 




8<S>{E V )S<P{E' V )) ^ 

Then the error of the extracted flux is obtained by integrating over a as 

((5*r M ) 2 ) s = j da((8^) 2 ) s P[a\D,I]. 

3. Model of inclusive neutrino-nucleus cross section and neutrino flux 


/J 2 dE v dE' v . (2.10) 


( 2 . 11 ) 


To examine the method outlined in the previous section, we construct pseudo data of neutrino- 
nucleus reaction by using a model of neutrino-nucleus reaction and neutrino flux. The neutrino flux 
is obtained by applying the formula of two-flavor oscillation as 


&(E V ) = 1—sin" 2023 sin" 


2-m „;„2 (Am.23L\ \ o 


V Ev ) 


W?v) 


(3.1) 


Here we used Am 2 3 = 2.5 x 10 3 eV 2 , sin 2 2023 = 10 and L = 295km. The initial flux flf( E v ) 


is 


taken from near detector neutrino flux of T2K experiment Ref. [10]. The solid(black) curve in Fig. 



Figure 1: Neutrino flux. The solid(black) and the dashed(red) lines show the neutrino flux <J> V and <f>y. 


1 is the neutrino flux we use to calculate the pseudo data, while the original flux is shown in 
dashed(red) curve for comparison. 

As a observable to test the method of flux extraction, we use double differential cross section 
of muon of neutrino induced inclusive reaction on 12 C, v u + 12 C —>■ /I ~ +X 


d 2 c 

dE^dO-n 


(3.2) 


In the energy region of neutrino flux in Fig. 1, the main reaction mechanism is quasi elastic scat¬ 
tering and pion production through the A(1232) excitation. The quasi elastic nucleon knockout 
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process is calculated using modified Fermi-gass model which incorporates the spectral function 
[12, 13]. For the delta excitation region we use the neutrino pion production amplitudes calculated 
in a formalism in which the resonance contributions and the background amplitudes are treated 
on the same footing [|l4|]. Detail on the reaction model is described in Ref. JTH]. The model of 
nuclear reaction is rather simple without taking into account the final state interaction, however the 
reliability of our formalism is tested against the electron nucleus scattering cross sections in the 
same theoretical framework; the calculated cross sections agree reasonably well with the existing 
data. Therefore, we believe it is good enough for the purpose of testing the method of neutrino flux 
extraction. 

We now are able to calculate the flux averaged double differential cross section cl 2 a / dE^clEl^ 
from the model of neutrino-nucleus reaction and neutrino flux <t> v in Eq. (3.1). 

d 2 c \ f d 2 c 


dE^ (-1 —l u 


= / dE v 


dE^dO.^ 


ME V 


(3.3) 


The pseudo data of double differential cross section is obtained as following. We assume muon 
energy resolution is lOOMeV and evaluate the average of cross section as, 

d 2 a 


d 2 o 

dE n dji 


y~i> 


i 

'e{^e[ 


d E^cl El n 


(3.4) 


We assume Gaussian distribution of the data with the mean given by Eq. (3.4) and the standard 
deviation by the 10% of the value of flux averaged double differential cross section at E u = I GcV. 
Fig. 2 shows the obtained pseudo data of the muon energy distribution at the muon scattering angle 

0,4 = 10 °. 



Figure 2: Pseudo data of neutrino induced inclusive muon production cross section on 12 C. The scattering 
angle of the muon is 0^ = 10°. 


4. Results and discussion 

We extract neutrino flux from the 19 points of pseudo data of the muon energy distribution 
shown in Fig. 2. They give D/ and (7/ of Eq. (2.5). D\ is calculated from the theoretical cross 
section of Eq. (3.2). 
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At first, we study the effects of the default model I on the extracted neutrino flux. We have 
examined two default models. The one is constant flux and the other is initial flux shown in the 
dashed(red) and solid(black) curve in the left panel of Fig. 3. The extracted flux using each default 



Figure 3: Default model dependence of extracted neutrino flux. The solid(black) and dashed(red) curves in 
the left panel show <f>® and constant flux as default model. The extracted flux corresponding to each 

default model is shown in the right panel. 

model is shown in the right panel of Fig. 3. One can see that, in this case, the extracted flux is very 
stable against the choice of the default model. 

We now compare the extracted flux d>J, // ;W (solid curve) with the expected ’exact’ neutrino 
flux <E> v (dashed curve) in Fig. 4. Extracted flux agrees well with the ’exact’ flux over the wide 
muon energy region. The extraction of the flux is only valid on the energy above the reaction 
threshold and the flux extraction did not work well near the threshold energy. We noticed that the 
extracted flux has structures in the 0.7 < E v < 1 GeV region. One can judge the structure is real 



Figure 4: Comparison of extracted flux <fdf £M (solid curve) and ’true’ flux <f> v (dashed curve). 

or not by estimating the error of extracted flux. In Table 1, the neutrino flux averaged in the peak 
region (0.28 < E v < 0.42 GeV), dip region (0.52 < E v < 0.70 GeV) and higher energy regions 
(0.72 < E v < 0.82 GeV), (0.86 < E v < 0.94 GeV) and (1.00 < E v < 1.08 GeV) are given. The 
second raw shows the calculated ’exact’ value of the average flux and the third raw shows the flux 
obtained from MEM together with the error calculated from Eq. (2.11). The central values of the 
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<pM EM agree well with the ’exact’ ones within the errors in all the energy regions. The size of the 
oscillation of seen in 0.7 < E v < 1 GeV is smaller than the estimated errors in that region. 

Therefore one can conclude that we sould not take the oscillation of the obtained flux seriously 
from the error estimation even without knowing the ’exact’ flux. 


£V(GeV) 0.28-0.42 0.52-0.70 0.72-0.82 0.86-0.94 1.00-1.08 

Exact 53.5 3.9 10.5 10.9 8.7 

MEM 56.5 ±6.5 5.1T1.7 14.0T3.7 11.6T3.7 10.7T1.9 


Table 1: Energy averaged neutrino flux. 

In summary, we have examined MEM to extract the neutrino flux from the inclusive muon 
production cross section. The results show the proposed method can be useful alternative tool to 
extract the neutrino flux provided accurate theoretical cross section is known. Further work to 
find the optimal observables for the efficient extraction of flux with the refined model of neutrino- 
nucleus reaction should be done in near future. 

The authors would like thank to Drs. M. Kitazawa and S. Nakamura for illuminating discus¬ 
sions. This work was supported by JSPS KAKENHI Grant Nos. 24540273 and 25105010. 
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